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We consider a spin coherent states description of a general quantum spin system. 
It is shown that it is possible to use the spin- 1/2 representation to study the 
general spin- J case. We identify the 1/2 spinor components as the homogeneous 
coordinates of the projective space associated to the complex variable that labels 
the coherent states and establish a relation between the two-component spinor and 
the bosonic Schwinger representation of a spin operator. We rewrite the equations 
of motion, obtained from the path integral for the evolution operator or partition 
function, in terms of the 1/2 spinor and define the effective Hamiltonian of its 
evolution. 



1 Introduction 

Coherent states provide a convenient way to represent quantum operators in 
terms of c-number functions. They have been used extensively in the liter- 
ature in various fields like quantum optics, condensed matter and quantum 
field theoryEl. They are closest to a classical description in the sense that they 
minimize the Heisenberg uncertainty relations. In-particular, spin coherent 
states have proved useful to deal with spin systems □. Since the algebra of the 
spin operators is more involved than that for boson or fermion operators, the 
construction of the spin coherent states is somewhat more complicated. Re- 
cently, we have pursued the similarities to the boson coherent states and have 
shown that it is possible to establish a closer parallelelism than previously 
realized. In particular, we obtained generalized representatives for spin opera- 
torsLJtl extending results obtained previously for bosonsu. Also, we considered 
the path integral for the transition amplitude and the partition functional and 
considered quantum Monte Carlo algorithms using coherent statesQ. 

The handling of spin operators is traditionally a difficult problem implying, 
in general, unphysical states thatJaave to be projected out. There are methods 
which do not include extra states™ but they are not easily generalized to spin 
values greater than one. Another advantage of the spin coherent states is that 
the spin value J enters only as a parameter and therefore all spin values can 
be handled similarly. The price to pay is a description in terms of continuous 
variables instead of a discrete number of states. The spin coherent states are 
therefore a convenient basis to perform a semi-classical expansion for large 
values of the spinllS .and, they appear naturally in the quantization of a semi- 
classical spin theorylIinI3. 
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The fact that J is only a parameter in the formalism suggests that one 
can reformulate the spin- J case in terms of the J — 1/2 case. In this paper we 
emphasize that indeed most results obtained for a general J-value are simply 
related to the same expression for J = 1/2 which is in general simpler to 
obtain. 

2 Spinor versus stereographic projection descriptions 

The spin coherent states can be defined by! 

\a>=R(6,<p)\JJ>= ) e aJ -\JJ>, (1) 
(1 + \a\ 2 ) J 

where R{0, (p) is a rotation through an angle 9 about the axis n — (sin ip, — cos cp, 0), 
normal to the z-axis and to the vector to — (sin 9 cos cp, sin 9 sin ip, cos 9) , and 
the state |JJ > is the top extreme eigenstate of J 2 and J z . The variable 
a = tan |e lv defines the stereographic projection from the south pole of the 
sphere to the plane passing through the equator, with complex coordinates 
a, a*. Under this projection, distances in the sphere are mapped into the 
plane according to (ds) 2 — d9 2 + sin 2 9d(p 2 — Ad\a\ 2 /(1 + |a| 2 ) 2 , and areas 
or spherical angles according to dil = sin 9d9dip = Ad 2 a/{1 + |a| 2 ) 2 , where 
d 2 a = ddlad^sa. Under the rotation considered above the vectors of the spher- 
ical basis are rotated according to (0, 0, 1) — * rh and 1/2(1, ±i, 0) — * h T , where 
2(1 + a*a)h- = (1 — a 2 ,i(l + a 2 ), —2a) and h + = h*_. These vectors satisfy 
to 2 = 1, h- ■ h+ = |, h± = and fh ■ h± = 0. The coherent states are not 
orthogonal and their overlap is 

(l + a'*a) 2J 

1 (i + ki 2 k(i + h 2 )'' (Z) 



^ (i + cose) " 1 J 



(3) 

where $ is the area of the spherical triangle with vertices to, to' and e z and 
is the angle between the two directions defined by to and to'. The coherent 
state \a > depends on the variables a and a*, which should be considered as 
independent. It is convenient to define matrix elements depending only on two 
complex variables (and not four) u. The same happens, in the case of bosons 
where the so-called holomorphic representation is usedllaEil. The easiest way 
to use this type of representation is to define new statesH ||a >= e aJ -\JJ > 
which depend only on a. Their overlap is < /3||l||a >= (1 + (3*a) 2J and the 
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decomposition of the identity becomes 



2J + 1 f 2 \\a >< a 



d a (ITww = L < 4 > 

The coherent states ||a > and < a|| for spin \ are given by the spinors 

A(a) = (') (5) 

A+(a) = (la*), (6) 

respectively, and one can write 

A t (a)A(a) = l + a*a (7) 
< a||i||a > = [At(a)A(a)] 2J . (8) 

The matrix element of an operator F is defined as F((3*, a) —< f3\\F\\a >. 
It depends only on a and /3*, and it does not depend on a* or (3. One can then 
use a* instead of /?*, and treat a and a* as independent variables, without any 
loss of generality. For some operators like the density matrix it is convenient 
to use the so-called diagonal representative given by the weight function of t 
decomposition of the operator in a superposition of coherent states projectorsl 
In connection with the holomorphic representation it is convenient to define 
the diagonal representative /(a, a*) given bya 

o 7 -l i r 

F=—— d 2 a\\a>f{a,a*)<a\\. (9) 



Since the spin coherent states are not eigenvectors of the angular momen- 
tum operators, it is important to know their action on them. These relations 



were 



obtained beforeQ 



e pj ~\\a > = \\at + p > 

e e' S +\\a> = (1+ P*a) 2J \\—^— > 

1 + p*a 

eP J >\\a> = e pj \\e- p a> . (10) 

The action of the infinitesimal operators on the coherent states follows imme- 
diatelyu and can be summarized as 

d 

J\\a >= {Jm + h-[(l + a * a )g^ ~ 2J&*}}\\a > (11) 
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The matrix elements of the operators J± and J z are given byuncHj < a\ J\a >— 
Jrn. Similarly, the diagonal representatives are given by j(a, a*) = (J + 
l)m. From eq. (11) follows that the coherent states minimize the uncertainty 
relation for any two spin components since, from the Schwarz inequality, in 
order that a state \<j> > minimizes the uncertainty relation, the states (J— < 
J, >)\4> >, with i = ±, 0, must be linearly dependent. 

The operators D(g) of the spin- J representation of the rotation group can 
be parameterized in different manners. The most intuitive one is to define the 
axis of rotation, defined by the unit vector rh and the angle of rotation 6. We 

have then D(g) = e l ^' J , where 6 = <j)fh. A specially useful set of parameters 
can be found looking at the form of these operators for J = 1/2, where they 
reduce to 



V{g) = e 



cos — V i— ■ asm — 

2 6 2 



fi —v 



(12) 



where /i = cos ^ + i cos 9 sin ^ and v — i sin 9e lip sin ^ . These two complex 
numbers are not independent, since they satisfy |/i| 2 + |^| 2 = 1. The inverse of 
an operator is obtained replacing fi by /x* and v by —v. Using these parameters, 
the disentangling theorem EJ is written as 

D(g) = e -£- / + e - 2In "* i *e7 Fi - (13) 
= et J -e 21n ^e~ir J + (14) 
The action of a general rotation on a coherent state is 

b{g)\\a>=^~v*af J \\ V -0^> . (15) 

The matrix element of the operator D{g) is then given by 

< a'||£>(s)||o! >= [a* + a'*v - v*a + n*a'*a] 2J (16) 

The matrix element of the spin 1/2 operator T>(g) between these spin | coher- 
ent states is given by 

T>(g; a* ,a) = fi + oi*v — v* a + /i*a*a 

= A\a)V(g)A(a). (17) 
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One concludes then that the matrix element for general spin J is simply 



<a\\D(g)\\a>=p(g;a*,a)Y J (18) 

i.e. the power 2 J of the corresponding spin i matrix element, with eq. (8) be- 
ing a special case. These results are easily understood if one uses the tensorial 

product of 2J spins i, J a , a — 1, . . . ,2 J, and write \ JJ >= \hh > ■■■ \~ > 

and (f.J = Y^ a =\4>-Ja- The matrix elements of a general rotation between 
coherent states factorize then in the product of the factors for each spin i, 
where the spin i coherent states are given by eqs. (5,6). 
The diagonal representative of D(g) is 

d(g; a, a*) = — — ■ — ^— 12f7+ir (19) 

which can be rewritten as 

J( g ;a,a*)= [p(g _ 1; ^ a)]2(J+1) . (20) 

i.e. the diagonal representative of D{g), for general spin J, is given by the 
power — 2(J+ 1) of the matrix element, for spin i, of the operator D~ 1 {g). 

A direct way to calculate the matrix elements of products of powers of 
spin operators consists in differentiating the appropriate generating function 
< a \\e uJa e vJb ■ ■■ \\a > and using the reduction of this matrix element to the 
2J power of the corresponding spin h matrix element. Repeated use of the 
identity <j a &b = S a b + ieabcO'c, where a a are the Pauli matrices, reduces the 
calculation to m a =< a\\a a \\a > / < a||l||a >. One concludes then that 
any correlation function is a function only of the vector J = Jfh. The same 
applies to diagonal representatives. In practice, however, correlation functions 
and cumulants, in particular, are most easily evaluated using the fact that the 
insertion at right (or at left) of additional spin operators into matrix elements or 
diagonal representatives of some operator is given by the differential operators 
Jm + /i_(l + a*a)-f^ and (J + l)m — h+(l + a*a)-^ (or their comfllex 
conjugates), acting on those matrix elements or diagonal representativesEj. 

If we consider the action of a rotation T>(g) on a general | spinor 

Ax(a)=(t) (21) 



we find that the relation between the components of the rotated A.' x (a) 
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f>{g)k\(a) and original spinor A\(a) can be written as 



a' 
A 7 



^ + M x 



(22) 



From eq.(15) it follows that the variable a, labeling the spin coherent states, 
has the same transformation law as the ratio a/A, of the two components of the 
spinor. This is the relation, found in geometry, between an affine space and 
the projective space associated with it. More precisely, we identify the 1/2 
spinor components as the homogeneous coordinates of the projective space 
associated to the complex variable which labels the coherent states. This 
coordinate transforms according to a bilinear or Mobius transformation. One 
recovers the well known isomorphism between the SU (2) group and the group 
of the bilinear transformations. The variable A is the rescaling necessary to 
have the first component of the spinor equal to one. 

Under a transformation in which the spinors A\(a) and A^(a) transform 
according to A\(a') = UA\(a) and A^(a') = A\(a)V, where U and V are not 
necessarily related, the Jacobian for the change of variables is given by 

(|A'| 2 ) 2 d 2 T7 = dct(7dct V(|A| 2 ) 2 d 2 ^. (23) 
A A 

The two components of the spinors transform separately as the numerator 
and the denominator of the bilinear transformations. This is very useful in 
practical calculations. For example, in the calculation of the integrals involving 
matrix elements or diagonal representatives, one can replace matrix elements 
between coherent states by matrix elements between spinors, and use rotations 
and not bilinear transformations, to simplify their evaluation. This, combined 
with the relation between the spin J and spin 1/2 quantities, allows a great 
simplification of the calculations. 

Also, the connection between rotations of vectors and rotations of spinors 
expressing the fact that the spin is a vector operator and given by TZ(rn) ■ a = 
Dm ■ aV^ 1 can be immediately obtained using the identity 

x - 2AAt * 
m(a) ■ a = -j-^ - 1 (24) 

in terms of the spinors A and At. Due to the homogeneity of the expression 
for m the vector J = J ^A\(rA\j / ^A^Aa^ coincides with the vector J = J fa 
previously defined, when we make the replacement a — » a/A, a* — > a* /\*. 



6 



3 Path integral and equations of motion 

The transition amplitude 

C7(a},t / ;a J ,t l )=<a/||T t e"^ j * l * mW \\a i > (25) 
can be obtained by the path integral expression i 

U(a},t f ;ai,ti) = J V 2 fi(l + a* f a(t f )) J e i ^ MC (1 + a*(t l )a l ) J (26) 
where 

h J , da* «.da , , « , . 
C = -il^-a^- a ^- H ^^ (27) 
is the Lagrangian, in which Ti. c (a* , a, t) is the normalized matrix element of 
the Hamiltonian. The dependence on a and a* is through m only, as explained 
above. The two boundary terms appear naturally from the construction of the 
path integrals and are necessary to have the boundary conditions a(ti) = ai 
and a*(tf) — a*j, when varying the action. Note that a(tf) is not the complex 
conjugate of a*c and that a*{ti) is not the complex conjugate of o^. They are 
the values taken by the variables a(t) and a*(t), satisfying these boundary 
conditions, at the times t = tf and t = ti, respectively. In particular, in the 
saddle-point approximation they are obtained solving the classical equations of 
motion for a(t) and a*(t) and evaluating their values at tf and ti, respectively. 
These terms have been shown to be necessary to obtain the exact solution using 
the saddle-point approximation in the case of a spin in a time-independent 
magnetic field, for instance a. In the discrete approximation the integration 
measure is V 2 fi = II^^ 1 2J ^~ 1 ( 1+ rf a "„ fc )2 ■ It can be seen as a consequence of the 

quantization of a classical problem with constraints 111. 
The usual transition amplitude is therefore given by 

< 0/ |T t e-* >= [V^ ( 1 + "f a MY e *X** fl + «-&M J 
il ' J ^\ l + a* f a f J \ 1 + ofOi J 

where af and a* are the complex conjugates of a*? and a^, respectively. The 
path integral representation for the diagonal representative of the evolution 
operator is similar 

It can be checked that, both for matrix elements and diagonal represen- 
tatives, the classical equations of motion, obtained varying the Lagrangian 
term and the boundary condition terms, are the classical limit o£the Dyson- 
Schwinger equations, with the appropriate boundary conditions 19. Also, the 



7 



Poisson brackets, or more precisely, the Dirac brackets for the quantization of 
a theory with constraints, of the matrix elements or diagonal representatives 
of any two spin components satisfy the commutator algebra of the quantum 
spin operators. 

The path integral expression for the partition function can be easily ob- 
tained and one finds using matrix elements 

Z = J V ^(t)e-*C dT[ - n ^ { ^ a - a *^ }+n ^> a)] (29) 

with periodic boundary conditions in the imaginary time. The expression for 
the case of the diagonal representative is similar u. 

Geometrically, the free term in the Lagrangian is the area in the sphere 
defined by the trajectory under consideration and by two great arcs from the 
north pole to the initial and final points of the trajectory. It can be rewritten 
as 

£ " = 4r x & ( 3 °) 

at 

where w = — %=r- , with u — el, reflecting the choice of the z axis as the 
quantization axis in the construction of the coherent states. In general, u can 
be any arbitrary unit vector. Changing it, the free Lagrangian is modified 

by SCq = 4t ( — x u- Su) i.e., the total derivative of a function of the 

initial and final positions (the difference in the areas defined by each of these 
points and the old and new u). Conversely, under a rotation, i.e., under a 
change of variables of the form a' — , the free Lagrangian changes by 

Cq — Co = j[i ln ^J"^,^ ) which is also the total derivative of a function of the 

initial and final positions. Finally, under a spin inversion, in which J — > — J, 

a — * — ^r, a* — * , the free Lagrangian transforms according to Co — > 

— C + ^(— iJln(^)) = — C — 2 J^jf , with the geometrical interpretation that 
the sum of the areas defined by a closed trajectory encircling the poles and 
viewed from them is 4ir, the total area of the sphere. 

The Lagrangian is totally rewritten in terms of the vector J. The integra- 
tion on the sphere can also be written as 1/ J 2 d? J5(J 2 — J 2 ). We arrive then 
at the path integral for spins 

Jvf6(f 2 -J 2 )e l I dtC (31) 

Similarly, a formulation involving only the spinors can also be found. Mak- 
ing the rescaling a — > a/A and a* — > a* /A* , the integration measure becomes 



S 



dfl = \\\ 2 d 2 a / (K^ K) 2 . In this expression the values of A and A* are arbitrary 
but fixed. In order to integrate over d 2 X one introduces two delta functions 
<5 (3ff A — A r ) and 5(SA — Aj). Using the identity 



5(f(x,y))S(g(x,y)) 



d(f,g) 



d(x,y) 



= 5(x - x Q )5(y - y ) (32) 

x=x ,y=y 



where xo and yo are such that f(xa,ya) = g{%o,yo) = 0, one is able to give 
other forms to the path integral. In particular defining 

W v /l + «> i Vx; 

where r(A*, A) is an arbitrary but not constant function, one findaJmmediately 
the path integral expression found by Jevicki and Papanicolaou LH 

Z= f\\V 2 a i ]\5[A\t)A{t)-2J]8{4>{t)]{4>,N} PB e i I dtC (34) 

derived using the Schwinger representation for spins and using the Dirac- 
Faddeev theory of constrained systems, extended to the case of second class 
constraints 113, to make the restriction to the physical space A 1 A — N = 2 J 
and where <j) = a\ + a\ is the gauge fixing condition. In terms of the function 
r(A*,A) this restriction and the gauge condition become \r\ 2 — 1 and r* = r, 
respectively. 

The classical equations of motion can be obtained varying the Lagrangian 
term and the boundary condition terms. In the case of the path integral for 
the matrix element, we find 

2J da _ dU 
1 {l + a*a) 2 ~dt ~ ~du* 

2j da* _ an 

(l + a*a) 2 dt ~ da ( ' 

with a(U) = ai and a*(tf) — a*? as boundary conditions, as expected0. Using 

(1 + a* a)dm / 'da = 2h + and its complex conjugated, the equations of motion 
can be rewritten as 

2 da dU 2/i_ 

in- 



(l + a*a) 2 dt dJ l + a*a 

2 da* OH 2h+ 
-in- — = — =; • (3d) 
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where J = J fa. In these equations H e ^ = is the effective magnetic field 
acting on the spin. Making the rcscaling a — > j and a* — * jt- and defining 
the effective Hamiltonian JC — H e ff ■ acting on the spinors, one finds that 
these equations of motion are equivalent to 

ih^± = K,A\ (37) 

and its hermitian conjugate. Since K. is an hermitian operator, the norm of the 
spinors, i.e. A* A A = A* A + a*a 1 is conserved. The equation of motion of the 
moment J, i.e., of the expectation value of the Pauli matrices in the spinor, is 

dJ dH f 

— = — - x J 38) 
dt dJ 

which is the equation of precession in the effective field ^j. 

The equation of motion of the spinor can be formally solved by 

A A (a) =U(t,t )A° x (a) (39) 

where 

U(t,t )=T t e~*fto dtK (40) 

is the evolution operator corresponding to the effective Hamiltonian K . One 
can take advantage of the connection between vector rotations and spinor 
rotations by defining the operator in spinor space 

J = J-a (41) 

having the property J 2 = J 2 . The equation of precession becomes 

in^ = [t,j] (42) 



and can be solved by 



3 = Uj^lJ- 1 (43) 



The case of many spins can be easily treated, repeating this procedure for 
each spin. The path integral follows immediately, with the coupling between 
spins arising from the Hamiltonian H. The saddle point equations follow from 
there and one finds the mean field approximation in which each spin i precesses 
in the effective field = ^j- resulting, in particular, from the interaction 

with external fields and the other spins. One introduces Pauli matrices <7j 
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for each spin. The effective Hamiltonian for each spin is ti = H- ff • f The 
total effective Hamiltonian is K. = YliJ^i- ^ is the sum of single spin terms, 
reflecting the character of the mean field approximation in which each spin 
precesses in the mean produced by the others. We note that in the treatment 
of continuous spin chains the matrix J is the central quantity forj-the Lax 
representation used by Takhtajanlia and Jevicki and Papanicolaou li^J in the 
derivation of nontrivial solutions of the soliton or instanton type. 

4 Summary 

We have replaced a description of a quantum spin operator in terms of classical 
variables (defined by the stereographic projection of the sphere from its south 
pole to the equatorial plane) by a spinor description (two-component spinor) 
and replaced the action of a rotation on the spin coherent states (which involves 
bilinear transformations of the complex variables) by the action of a rotation 
(matrices operating in spinor space) on the two-component spinor. This refor- 
mulation appears naturally for J = 1/2 but we have shown that it generalizes 
in a straightforward way to all spin values. We have noted that this spinor 
description is related to the description of a spin system by Schwinger bosons. 
Both descriptions enlarge the space and a projection to the proper subspace 
is required. We have reexpressed the equations of motion in terms of the 1/2 
spinor and obtained simple expressions that have a formal appearance familiar 
from standard dynamics of two-level systems. This enables a simpler treat- 
ment of a general spin- J system. Interacting spin systems are dealt with more 
conveniently in terms of an evolution operator that acts on spinor space and 
in terms of a transfer matrix that connects different points in real space. Also, 
we have made contact to previous treatments where the Lax representation is 
used in the context of the search for nontrivial solutions of the saddle-point 
equations via the inverse scattering method. 
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